Electronic and optical spectra in a diluted magnetic semiconductor multilayer 
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The effects of random distribution of magnetic impurities with concentration x in a semiconductor 
alloy multilayer at a paramagnetic temperature are investigated by means of coherent potential ap- 
proximation and tight-binding model. The change in the electronic states and the optical absorption 
spectrum with x is calculated for weak and strong exchange interactions between carrier spins and 
localized spin moments on magnetic ions. We find that the density of states and optical absorption 
are strongly layer-dependent due to the quantum size effects. The electronic and optical spectra 
are broadened due to the spin fluctuations of magnetic ions and in the case of strong exchange 
interaction, an energy gap appears in both spectra. Furthermore, the interior layers show higher 
contribution in the optical absorption of the system. The results can be helpful for magneto-optical 
devices at a paramagnetic temperature. 



I. INTRODUCTION 



Semiconductor multilayers are materials with attrac- 
tive electronic and optical properties and possible appli- 
cation in nanodevices [l|. Such unique properties which 
are caused by their dependence on the number of atoms 
in a confined direction, have led to a rapid growth on 
research activities in this field in recent years. For in- 
stance, it has been shown that for semiconducting multi- 
layer armchair graphene nanoribbon with more than two 
layers, semiconductor-to-metallic transition is obtained 
with increasing electric field, which can be utilized to en- 
hance the eflticiency of graphene based transistor devices 
0. In addition, for the thicker systems (i,e. more lay- 
ers) , this transition occurs at a smaller electric field which 
makes them to enhance the on/off ratio of a field-effect 
transistor device [1]. 

Furthermore, magnetic and nonmagnetic impurities in 
layered structures have received considerable attention 
[3l-lll|. Spin-dependent scattering of charge carriers in 
magnetic multilayers is the origin of magnetotransport 
effects and has been the subject of intensive interest in 
recent semiconductor electronics research [i-|il|. Doping 
of magnetic ions in nonmagnetic semiconductors (NMSs) 
introduces not only magnetic moment but also free car- 
ries in the system. In diluted magnetic semiconductors 
(DMSs), the host material is an alloy semiconductor of 
the type II- VI or III-V semiconductor with characteristic 
of AB, in which the atom A is substituted by a magnetic 
atom with concentration of x [H, [l3|- For example, Mn 
ions as magnetic impurities into GaAs and InAs act as ac- 
ceptors in the III-V compound semiconductors [H, [111 . 
These semiconducting alloys are a novel class of ferro- 
magnetic materials and one of the promising materials for 
spintronics [Tsl . The exchange interaction between a car- 
rier spin and localized spin of magnetic ions, plays an im- 
portant role in magneto-optical effects in DMSs jl3, 13 ■ 
Because this interaction strongly affects on band splitting 
that is observed in magneto-optical measurements such 
as magneto-reflectivity and magneto-absorption spectra. 



Thus, the most powerful tool for studying the exchange 
interaction between a carrier and localized spins in DMS- 
based structures is optical measurement. Furthermore, 
because of the ternary nature of DMSs, one can tune 
the lattice constants and band parameters by varying the 
composition of the material. Therefore, most of the mag- 
netic semiconducting alloys are excellent candidates for 
preparation of quantum wells, superlattices, and other 
heterostructures that involve band gap engineering 

In previous studies [2ll - l23l | , by applying the single-site 
coherent potential approximation (CPA) [2J| to semicon- 
ductor alloy multilayers, the effects of chemical (non- 
magnetic) disorder on optical absorption spectrum were 
investigated. The results showed that, the photon ab- 
sorption by the system is layer dependent and strongly 
changes by variation of the scattering-strength of the 
chemical impurities 21|. In the present study, we ex- 



tend our previous approach [2l| for the magnetic impuri- 
ties at a paramagnetic temperature, and investigate the 
electronic and optical properties of DMS multilayers. To 
study the optical absorption spectrum in the presence of 
the magnetic impurities, we use the Onodera-Toyozawa 
theory [255, which is related to a binary mixed crystal or 
an alloy consisting of atoms (or molecules) A and B, in 
which the transition dipole moments of A and B atoms 
are equal. We should point out that, the application of 
the Onodera-Toyozawa theory together with the CPA for 
DMSs was first employed by Takahashi [26] in which the 
optical band edge was studied for bulk materials. The 
calculations, however, were performed only for a simple 
semicircular model density of states. 

In this study, owing to the substitutional disorder of 
the magnetic ions in a DMS multilayer, the effect of car- 
rier scattering due to the presence of both magnetic and 
nonmagnetic atoms will be studied within the single-site 
CPA. Due to the quantum size effect in the multilayer 
systems the electronic states at each layer strongly de- 
pend on the location of that layer relative to the sur- 
face layer. Hence, each layer has a different response to 
the carrier scattering process in the layer and should be 
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taken into account in the CPA equations. Since in a low- 
dimensional system many important physical properties 
and characteristics such as optical properties, are gov- 
erned strongly by the electronic structure of the system, 
then the local density of states (LDOS) will be analyzed 
precisely. Therefore, the dependence of electronic and 
optical spectra on the carrier energy, magnetic impurity 
concentration, and the strength of exchange interaction 
between the carrier spin and the spin of magnetic ions 
will be studied with more details. This paper is orga- 
nized as follows: we start with introducing the model 
and formalism in Sec. II. In Sec. Ill, the results of the 
numerical calculations for DMS multilayer are given. The 
conclusion is presented in Sec. IV. 



II. THE MODEL 

In this study, we consider a semiconductor multi- 
layer described by the single-band tight-binding model 
with the nearest-neighbor hopping and the on-site delta- 
function-like potential on a simple cubic lattice in which 
one of the dimensions (the z direction) is confined. The 
number of layers along the confined direction is and 
the label of each layer is denoted by n, hence 1 < n < A^^. 
The multilayer is a semiconducting alloy of the form 
Ai-xMxB, where the sites of A atoms such as Ga can 
be occupied by M atoms such as Mn with concentration 
X and B atoms such as As remain unchanged. Due to 
the low carrier density in the system, the interaction be- 
tween carriers is ignored in the calculations. Accordingly, 
the one-electron or a Frenkel exciton Hamiltonian in this 
system is given by: 



H 



2^2^2^trn.r'n'\r,n,ti){r',n',fi\ (1) 



r,n r',n' 



where, |r, n, fi) is an atomic orbital with spin /i(=t or 
at site (r, n). Here, r denotes the positional vector 
in the x — y plane of the layer n. The hopping integral 
trn.r'n' IS equal to t for the nearest neighbor sites and 
zero otherwise. We assume that the hopping integral 
only depends on the relative position of the lattice sites. 
Therefore, the type of disorder in our model is considered 
to be diagonal and accordingly, the off-diagonal disorder 
is outside of the scope of the present study. The random 
site energy u^i^ is assumed to be and with 
probabilities 1 — a; and x when the lattice site (r, n) 
occupied by the A and M atoms, respectively [21I.I26 
For A and M sites we have: 



£A^\r,n,fi){r,n,fj,\ 



(2) 



,M 



k, n, H) [£m5^u - ISrn ' CT^i/] (i", n, v\ (3) 

where JSm • tr^i/ is the fc-independent exchange interac- 
tion in which | Sm | = S* is the local spin operator of the 



M atom and a is the Pauli matrix for carrier spin. In 
this study, instead of taking the quantum fluctuation of 
the localized spin of the magnetic ions into account, the 
classical spin approximation (i.e. S — >■ oo) is used, while 
IS is a constant parameter [2^, [23| • 

According to the CPA, a disordered alloy can be re- 
placed by an effective periodic medium in which the po- 
tential of all sites is replaced by an energy-dependent 
coherent pot ential, except one site which is denoted by 



impurity |24{ . Therefore, the multilayer Hamiltonian H 



is replaced by the following effective medium Hamilto- 



nian: 

HcS = EEE[Sn(w)(5r,r"5ri,ri' - i] |r, n, ^) (r', n', ^| 

rn r'n' fJ. 

(4) 

where S„(a;) is an energy-dependent self-energy and 
called the coherent potential. Because of the absence 
of translational symmetry of the system along the z- 
direction, this self-energy depends on the layer number n. 
Furthermore, the self-energy is spin independent, because 
the system is at a paramagnetic temperature (T 3> Tc). 

The physical properties of the real system can be ob- 
tained from the configurationally averaged Green's func- 
tion, (G')av = ((^^ — -ff)^^)av which is replaced in the 
CPA with an effective medium Green's function G — 
(u! — Hcff)^^ ■ A direct consequence of this replacement 
is the fact that the effective scatte ring of a carrier at the 
impurity site is zero, on average j24| . The matrix ele- 
ments of the effective Green's function are calculated by 
the following Dyson equation [23| : 



G„,„,(r,r') = G"„,(r,r') + 



N, 
n" — 1 r" 



n ^ n 'n 



(r",r') 



(5) 



where the variable uj has been suppressed for simplicity. 
Here, G^{uj) is the clean system Green's function and its 
matrix element is given by 27j : 



^2 

G°,„'(r,r') ^ 



e=i 



IBZ 



(6) 

In this equation, £ denotes the subband number (mode), 
k|| is a wave vector parallel to the layer, and the in- 
tegral is taken over all the wave vectors in the first 
Brillouin zone (IBZ) of the two-dimensional lattice [2ll |. 

„, {£, k|| ) is the mixed Bloch-Wannier representation 
of G" which can be expressed as 



G" '(^,k||;c.) 



hn,n' ('^) 



UJ + ir] - e£(k||) 



(7) 



where /i„,„'(^) = (jv^+i) sin(jv^n) sm{j^n'), 77 is a 
positive infinitesimal and ef(k||) is the electronic band 
structure for mode £ (211] . 
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To obtain the CPA condition, we define a perturbation 
potential energy as 



(8) 



where Vr^n = ^ ^it^ ^^"^ ^r,n = ^^^r the 

M site are given by the following equations 27]: 



(9) 



I r,n, ^)[(eM -S„)(5^,,. - /Sr 



(10) 

Substitution of a fraction x of the element A by mag- 
netic impurity atom M will include both substitutional 
disorder and spin scattering, hence we can expect a semi- 
conductor with magnetic properties. In such a case, the 
CPA equation for the system is expressed as: 



, (11) 



where (tr'n) represents the complete scattering asso- 
ciated with the isolated potential v^^^ {v^^n) due to the 
A{M) atom in the effective medium and (• • • )spm denotes 



average over the spin scattering at the M site [27[. At 
the paramagnetic temperature the orientation of local- 
ized spin is completely random; hence the probability of 
each state is 1/2. Therefore, the associated t-matrices for 
an arbitrary site r in the nth atomic layer are given by: 



EA - S„(a;) 



and 



,M 



1 - {eA~^n{io))Fniuj) 
£m - is - i;„(a;) 



i-{eM-is-i:,M)F^H 

Em + IS - T.n{uj) 



(12) 



(13) 



(£M + /^-S„(w))K(w) 



where F„{uj) = G„.„(r, r; oj) is the diagonal element of 
the effective Green's function matrix [^J, • From the 
CPA condition, Eq. ([TT|) . one can derive an equation 
for the self-energy of rtth layer. Then, using such an 
equation and also Eq. ([5]) which gives a system of lin- 
ear equations, one can obtain self-consistently the self- 
energy E„(a;) and the Green's function Fn{uj) in each 
layer. Then, the LDOS per site in the nth layer, gn{Lu), 
is calculated by 



gn{uj) = --ImF„(i:j). 

TT 



(14) 



On the other hand, to calculate the optical absorption 
spectrum, we assume that both the magnetic and non- 
magnetic atoms have equal transition dipole moments. 



Therefore, when the £th subband is optically excited, the 
layer- dependent optical absorption is given by the contri- 
bution of k|| = component of the electronic states that 
is related to the F point in the IBZ \2l\. Accordingly, the 
optical absorption of the nth layer, due to the creation 
of an exciton in the system, can be defined as 



1 

An{uj) = ImVG„,„(^,k|| =0;u) 

TT ' * 



(15) 



Note that, the calculation of A„(a;) is not restricted to 
the lowest £ = 1 subband, because the numerical results 
showed that the contribution of k|| = component in the 
other subbands, i.e. ^=2,3,4 and 5, is finite and hence, 
all the subband contributions must be included. In the 
next section we present the numerical results of LDOS 
and optical absorption spectrum for the system under 
consideration. 



III. RESULTS AND DISCUSSION 

Based on the above formalism we study the elec- 
tronic states and the optical absorption spectrum in a 
DMS alloy multilayer consisting of five atomic layers, i.e. 
Nz = 5. All the energies are measured in units of t 
and we set eA=0 as an origin of the energy. We per- 
form the calculations for two sets of weak {em — —0.75 1, 
IS = -1.5 1) and strong (em = -4.0 IS = -3.5 i) in- 
teractions. It has been shown that, the band gap opening 
in DMSs depends on the value of magnetic ion chemical 
(or spin independent) potential and the strength of ex- 
change interaction [23. Therefore, by choosing the above 
values for sm and IS we can model the system to show 
gap opening in the electronic states due to the doping 
of M atoms into the AB structure. Since the system 
consists of five atomic layers, one can expect five energy 
subbands in the band structure of the system [27|]. Each 
subband is attributed to one of the atomic layers. In ad- 
dition, due to the geometrical symmetry of the system 
in the z- direction, the electronic states and hence the 
optical properties in the layers n = 1 and 5 and also in 
the layers n = 2 and 4 are the same. Therefore, we only 
present the LDOS and the optical absorption spectrum 
for the layers n = 1, 2 and 3. 

Fig. 1 shows the dependence of electronic states on 
the carrier energy uj and the alloy concentration x in the 
case of weak interaction, i.e., Em = —0.75t, IS = —1.5t. 
Because of the 2D nature of the atomic layers of the sys- 
tem, one can see the van Hove singularities and steplike 
features in the LDOS which strongly depend on the mag- 
netic ion concentration. The energy position of these fea- 
tures is different in various layers due to the quantum size 
effects which cause different electronic states in different 
layers. In the case of a; = the system is in completely 
nonmagnetic case, i.e., the original AB system, while in 
the case of a; = 1, the AB system changes to the MB 



(a) 



(b) 



(c) 




(a) 



3^ 
<" 



(b) 



< 



(c) 



CO 

< 



FIG. 1: The LDOS for a DMS multilayer with em = -0.75 f 
and IS = — 1.5t as functions of energy u) and alloy concen- 
tration X for the atomic layers (a) n = 1, 5, (b) n = 2, 4, and 
(c) n = 3, respectively, oj is measured in units of t. 



system that is a clean magnetic system. When the al- 
loy concentration is increased, the bottom of the bands 
shifts towards lower energies while the top of the bands 
remains approximately at a fixed energy. The shift of the 
bands is due to the enhancement in the spin fluctuation 
of M ions with increasing x. In the limit of x — 1, a dip 
appears close to the center of the band that for x = 1, 
its energy position is at £ = — 0.75t, i.e. at the on-site 
energy of the M ions. The optical absorption spectrum, 
associated with the above electronic states, is shown in 




FIG. 2: The optical absorption spectrum for a DMS mul- 
tilayer with Em ~ — 0.75t and IS = —1.5 1 as functions of 
energy lo and alloy concentration x for the atomic layers (a) 
n — 1,5, (b) n = 2,4, and (c) n = 3, respectively, uu is 
measured in units of t. 



Fig. 2. The results show sharp features in the optical 
spectra corresponding to some of the features seen in the 
LDOS at low energies. These features are remarkable 
in the case of the clean system, i.e. in the AB system. 
When the magnetic ions are introduced, that is in the 
case of alloy system, the absorption peaks are broadened 
and the height of sharp features decreases. Because of 




FIG. 4: The same as Fig. 2 but for bm = -4.0 1 and IS = 
~3.5t. 

the dependence of A„(aj) to the electronic states through 
the imaginary part of Green's function, one can see a 

shift in the optical spectra toward lower energies, similar the A and M atomic states are amalgamated into a single 

to that of the electronic states. An important feature band when the exchange coupling is small, 

in Fig. 2 is the dependence of optical absorption spec- To study the effect of strong exchange interaction on 

trum on the layer number. This means that the photon the electronic and optical properties of the system, we 

absorption in low dimensional systems is different in var- have shown in Fig. 3 and 4 the layer-dependent of LDOS 

ious layers. In addition, one can find that in the case and optical absorption for e = —A.Qt and IS — —Z.bt. 

of weak interaction no energy gap appears in the optical The influence of magnetic impurity causes an impurity 

spectrum. Accordingly, the density of states and the ab- band at the bottom of the host band, whose bandwidth 

sorption spectrum in the weak interaction indicate that depends on the M concentration x. With increasing x. 



the contribution of spin scattering by the locahzed spins 
in the scattering process of the carriers increases and the 
impurity band and hence, the whole electronic spectrum 
is broadened. This broadening which does not appear 
in LDOS of the NMS alloys 2l| is a consequence of the 
localized spin fluctuation of magnetic impurities. It is 
evident that the strong exchange interaction splits the 
host band into two subbands. We see that, the behavior 
of the electronic states of each impurity band completely 
depends on the magnetic ion concentration. For x=l, the 
two bands are symmetric and correspond to the parallel 
coupling and antiparallel coupling between the carrier 
spin and the localized spin, in agreement with the results 
of the CPA for ordinary magnetic semiconductors at a 
paramagnetic temperature [28l. [29^ . The energy position 
of the electronic band center is located at the energy 
£ = —A.Ot, i.e. at the on-site energy of the M ions, 
which is a consequence of the classical spin treatment of 
the localized spins at a paramagnetic temperature. Note 
that, the quantum effect of the localized spins introduces 
some further complications in the results and the band 
symmetry is broken. 

One of the main differences between the electronic 
states in the magnetic and nonmagnetic semiconductor 
alloys is the variation of with increasing x. In the 

CPA method and in the absence of magnetic impurities 
[U, a NMS alloy for a; = 1 is converted to the same 
original NMS but with an energy shift in the electronic 
states. This means that, the behavior of LDOS with en- 
ergy variation at a; = 1 is the same behavior of the LDOS 
at X = 0. By doping a NMS with magnetic impurities, 
however, the behavior of LDOS might be completely dif- 
ferent with the electronic states of the original undoped 
system, specially at x = 1. This difference appears be- 
cause of the multiple scattering effect due to the exchange 
interaction between the carrier spin and the localized spin 
of magnetic ion, which is most remarkable in both cases 
of a; = and x = 1. Because of the dependence of op- 
tical spectrum on the electronic states, this remarkable 
difference for the cases of a: = and a: = 1 can be clearly 
seen in this study. It is important to note that, if we set 
IS— 0.0, the present formalism will give the same results 
for a NMS multilayer in which there is no band broaden- 
ing in the electronic spectra [21 1. 

The electronic band splitting in the strong exchange 
interaction regime creates the two subbands in the op- 
tical absorption spectrum. For < a; < 1, the lower 
energy band corresponds to the case that the carrier spin 
is parallel with the localized spin, while the higher band 
belongs to the A atom and the case in which the carrier 
spin is antiparallel with the localized spin. Therefore, 
such an optical absorption spectrum might correspond 
to the persistence type [1^. For a; = 1, the spectrum 
completely belongs to the M atoms, i.e. the clean mag- 
netic system. Optical band broadening in the all layers 
can be clearly seen in the spectrum due to the chemi- 
cal and magnetic disorder. Furthermore, the bands show 
a sharp feature at the optical band edge, depending on 




FIG. 5: The total optical absorption spectrum for the two 
DMS multilayers as a function of energy u and alloy concen- 
tration X. In (a) Em = — 0.75i and IS = —1.5t and in (b) 
Em = — 4.0t and IS = —3.5 1. to is measured in units of t. 



the value of x. The strength of this feature increases for 
the interior layers and indicates that most of the pho- 
ton energy can be absorbed by the interior layers of the 
system. Roughly, for x < 0.5, the sharp features ap- 
pears at the bottom of the higher band, while for x > 0.5 
the respective features appear at the bottom of the lower 
band. Note that, the band broadening, which is accom- 
panied by energy shift of the bottom of the optical bands, 
demonstrates the operation of a magnetic semiconductor 
multilayer as an optical device which can be tuned at 
desired wavelength. 

Another set of the results that we present here is the 
bulk (total) optical absorption for the weak and strong 
exchange interactions. Due to the fact that the optical 
absorption of each layer is not experimentally measurable 
yet, we have calculated the total magnitude of this quan- 
tity here as well (see Fig. 5). For this purpose, we make a 
summation over all five layers to obtain the total optical 
absorption of the system, i.e. A{uj) = '^n=i^ri{^)- All 
features that we see in Figs. 5(a) and (b) are compatible 
with those of the individual layers. It is clear that the 
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sharp features of the amalgamation type (Fig. 5(a)) at 
the bottom of the optical band are much stronger than 
those of the persistence type (Fig. 5(b)). The results also 
show that the optical absorption of a DMS multilayer has 
an obvious broadening in comparison with the NMS mul- 
tilayers [111 , and in the case of larger IS the amount of 
broadening is higher and a gap opening occurs. There- 
fore, the transition from the amalgamation type to the 
persistence type in the DMS multilayers depends on the 
strength of the exchange coupling IS, and the magnetic 
ion concentration x. 



IV. CONCLUSION 

We have studied the effects of magnetic impurity on the 
electronic and optical properties of a semiconductor mul- 
tilayer under the assumption that the transition dipolc 
moments of the magnetic and nonmagnetic ions are the 



same. Using the single-site CPA for a random distri- 
bution of the impurity atoms, we investigated the influ- 
ence of exchange interaction strength and the impurity 
concentration on the LDOS and the optical absorption 
spectra of the system. In such multilayers, the magnetic 
impurities shift the spectra towards lower energy-side, 
broaden the bands, and strongly affect on the sharp fea- 
tures of the spectra, even in the case of weak exchange 
coupling. In the case of strong exchange coupling, an en- 
ergy gap opens in the electronic and optical bands. The 
interior layers show a higher contribution in the optical 
absorption process of the system relative to the surface 
layer, together with the fact that the amalgamation type 
alloys in comparison with the persistence type alloys need 
a higher amount of photons energy to make an excita- 
tion. These results might be helpful for the development 
of DMS multilayers in magneto-optical devices at a para- 
magnetic temperature. 
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